o 
o 
o 

(N 






C/3 



o 



o 



X 



Collapsed and adsorbed states of a directed polymer chain in 

two-dimensions. 

Pramod K Mishra and Yashwant Singh 



Jj^ [ Department of Physics, Banaras Hindu University, 

CN ; Varanasi 221 005, India 

(N 



(February 1, 2008) 



Abstract 



C^ , A phase diagram for a surface interacting long flexible partially directed poly- 

a. 

mer chain in a two-dimensional poor solvent where the possibility of collapse 



in the bulk exists is determined using exact enumeration method. We used a 
model of self attracting self- avoiding walk and evaluated 30 steps in series. An 



1^ \ intermediate phase in between the desorbed collapsed and adsorbed expanded 

in 

■^ ' phases having the conformation of a surface attached globule is found. The four 

cn 

(^ 
o 
o . 

f^ I panded, {iv) surface attached globule are found to meet at a multicritical point. 



phases, viz. (i) desorbed expanded, (ii) desorbed collapsed, (Hi) adsorbed ex- 



C^ ! These features are in agreement with those of an isotropic (or non directed) 

polymer chain. 
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I. INTRODUCTION 

The subject of adsorption of a long flexible polymer chain immersed in a poor solvent on 
an impenetrable surface with an attractive short-range attraction has received considerable 
attention in recent years Kfj-Q. This is because an attractive surface may lead to adsorption- 
desorption transition from the state when the chain is mostly attached to the surface to the 
state of detachment when the temperature is increased. This behaviour finds applications in 
lubrication, adhesion, surface protection, etc. 

It is known that when a chain interacts with a surface its conformational properties get 
modified in comparison with its bulk properties 0,@]. This is because of a subtle competition 
between the gain of internal energy and the corresponding loss of entropy due to the constraint 
imposed on the chain by the impenetrable surface. This competition leads to the possibility 
of co-existence of different regimes and multicritical behaviour. 

The essential physics associated with the behaviour of a surface interacting polymer chain 
in a good solvent is derived from a model of self-avoiding walk (SAW) on a semi-infinite 
lattice. If the surface is attractive, it contributes an energy e^ (< 0) for each step of the 
walk on the surface. This leads to an increased probability characterized by the Boltzmann 
factor uj = exp(— e^/fc^T) of making a step along the surface, since for e^ < 0, cu > 1 for 
any finite temperature T {kp is the Boltzmann constant). Because of this the polymer chain 
get adsorbed at low temperatures on the surface while at high temperatures all polymer 
conformations have almost same weight and non adsorbed (or desorbed) behaviour prevails. 
The transition between these two regimes is marked by a critical adsorption temperature T^, 
with a desorbed phase for T > Ta and adsorbed phase for T < T^ as shown in Fig. 1. At 
T < Ta the chain acquires a conformation such that a fraction of monomers get attached to 
the surface and others form a layer of finite thickness parallel to the surface. The thickness of 
the layer diverges at T = T^. One may define the crossover exponent 0, at T = Ta as Ng ~ A^''^ 
where A^ is the total number of monomers and A^s the number of monomers on the surface. The 
transition point T^ is a tricritical point [H. Both the surface and the bulk critical exponents 



have been calculated using renormalization group methods ^ , transfer matrix methods 
exact enumeration method M and Monte Carlo methods B. For a two dimensional system 



exact values of the exponents have been found by using conformal invariance |T^ 




Fig. 1; Schematic representation of polymer conformations in (a) desorbed and (b) adsorbed 
states. 

The situation is, however, different when the surface interacting polymer is a in compact 
globule state. In the vicinity of an impenetrable attracting surface the monomer-monomer 
attraction which is responsible for the collapsed state of the chain and the surface-monomer 
attraction responsible for the adsorption will compete. This competition may give rise new 
features in the conformational behaviour of the chain. For instance, a globule can get attached 
to the surface at some value of monomer-surface attraction without any significant modifica- 
tion in its structure as shown in Fig. 2. Such a phase has been shown to exist in between 
the bulk collapsed and adsorbed phases in both two and three dimensions for isotropic i.e. 



non-directed) polymer chain and is called surface attached globule (SAG) phase [jTI],|T2|. In 
an adsorbed phase the number of monomers on the surface are proportional to A^ and the 
thickness of the layer formed by the chain parallel to the surface proportional to A^°. On the 



other hand the bulk globule has one correlation length (the radius of gyration) which scales 
with N as R (X N^/'^ where d is the dimensionality of the space. In the SAG phase we have 
two correlation lengths, one parallel and another perpendicular to the surface. They may scale 
with N as -Rii oc A^^n and R\_ oc A^^'-^ where ^ < a;ii < -rr and < a; i < ^. 




Fig. 2; Schematic representation of adsorption of a compact globule; (a) shows a desorbed 
globule, (b) and (c) surface attached globule and (d) the adsorbed state. 



In this article we investigate the surface adsorption and collapsed of a directed polymer 
chain on a square lattice. Since for two dimensional system surface is a line, it can have only 
one adsorbed phase which we refer to as adsorbed expanded (AE) phase. The problem of 
adsorption and collapse of directed polymer has been studied using transfer matrix technique 
||r3| . The phase boundary of the collapsed phase separating it from the adsorbed expanded 
and the bulk expanded has been obtained exactly |]14| as well as numerically [0. These cal- 



culations show that there are three phases, viz. desorbed expanded (DE), desorbed collapsed 
(DC) and adsorbed expanded (AE). These three phases meet at a multicritical point. In these 
calculations a condition xu = 1 (where x is the step fugacity and u the Boltzmann factor for 
the monomer-monomer attraction) is used. This condition is valid only at the collapsed phase 
boundary. Therefore only the phase boundaries separating the DC phase from the AE and 
DE phases were obtained by these calculations. The phase boundary separating the DE phase 
from the AE phase has been obtained numerically by Veal et al |l^ . Since the transition from 



the DC phase to the SAG phase takes place in the region where xu > 1 the calculations based 
on the transfer matrix could not show this transition. 

II. CALCULATIONAL DETAILS 

We consider partially directed self-attracting self-avoiding walks (PDSASAWs) on a square 
lattice restricted to quarter space. Walk starts from the corner of the surface. Let Cn,Ns,Np be 
the number of PDSAWs with A^ steps, having Ng (< A^) steps on the surface and Np nearest 
neighbours. We have obtained C^^Ns^Np for A^ < 30 for square lattice by exact enumeration 
method. We prefer this technique because in this case the scaling corrections are correctly 
taken into account by a suitable extrapolation scheme |^. 

We associate energy e^ with each monomer on the surface and e^ for the monomer-monomer 
interaction. Partition function of the attached chain is 

Zn{uj,u)= Y. Cm,n.,NpUj''^u''^ (2.1) 

Ns,Np 

where a; = e"'^^/'^^ and u = e"'^^/'^^. to > 1 and u > 1 for attractive force. Reduced free energy 
for the chain can be written as 

G{uj,u) = lim — log Ziy{u,u) (2.2) 

N^oo Jy 

In general it is appropriate to assume that as A^ — > oo 

Zn{uj, u) ~ A^^-V(^, u)^ (2.3) 



where ^{uj,u) is the effective coordination number and 7 is the universal configurational 
exponents for walks with one end attached to the surface. The value of fi{uj, u) can be 



estimated using ratio method |T6| with associated Neville table or any other method such as 



Pade analysis IT^ or differential approximants [|T^,0. 
From equations (2.2) and (2.3) we can write 

log/i(a;,M)= lim -— log Zjq^uj^u) = G{ui,u) (2.4) 

N^oo Jy 

Zn{u!,u) is calculated from the data of Cn,Ns,Np using equation (2.1) for a given u and 
u. From this we construct linear and quadratic extrapolants of the ratio of Z]\}{uj,u) for 
the adjacent values of A^ as well as the alternate one. Results for alternate A^ give better 
convergence. When u = 1 and u = 1 the value of /i is found to be 2.39. 

The surface critical exponent 71 can be calculated using the relation [Q 

^ ^'~ log(iV/iV - 2) ^^-^^ 

and the end-to-end distance exponent u using the relation 



l \og{Rl/Rj 

2 log(A^/A^ - 2) 



^^xxus,..^/..^_,; ^2.6) 



where the superscript "0" indicates that the quantities correspond to the bulk without the 
surface and R^ is the end-to-end distance of a chain of A^ monomers. 

The value of udu) at which polymer gets adsorbed/attached to the surface for a given 
value of u is found from the (i) plot of Gs{uj, u) = G{uj, u) — G{u, u = 1) which remains equal 
to zero until uj = uOc and increases consistently as a function of uo for uj > ujc, (ii) from the 
plot of d'^Gs{uj, u)/des'^ at constant u and {in) from the plot of 7° — 71 as a function of uj for 
different A^. The value of uJc found from the plot of Gs{uj, u) is slightly lower than the peak 
value of d'^Gs/des'^- It is, however, observed that as A^ is increased from 22 to 30 the peak 
value shifts to smaller u and appears to converge on the value of u found from Gs{uj, u) plot. 
We therefore choose the value of Uc found from the plot of Gs{uj^ u) and determine lines uJc{u) 
and ujci{u) (see Fig. 3) by this method. For m = 1, the value of uJc is = 1.70 which is in very 
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good agreement with the value 1.707 reported in ref. 0. In another approach equation (2.5) 
is used to calculate 7° — 71 for different A^ and the values are plotted as a function of u. The 
value of ujc is found from the intersection of successive approximation to 7° — 71 in the limit 
N —>■ 00. This method, however, fails for u values above the 6'-point, i. e. where the desorbed 
phase is in a collapsed globule state. Below the 6'-point both the methods yield same values 
of Uc which are in very good agreement. 

The phase boundary separating the expanded and collapsed phases is calculated from the 
plot of Gb{uj, u) = G{uj, u) — G{uj, M = 1) as a function of u for a given uj. Here G?, measures 
the energy arising due to monomer- monomer attractions. This is however, not zero in the 
expanded phase and, therefore, cannot be used as an order parameter for the expanded and 
collapsed transition. We, however, find that the value of Gb{uj,u) as a function of u shows 
sudden rise at certain value of u. We locate this point from the peak of d'^Gb/dtp which 
corresponds to a specific heat peak and shows divergence as A^ -^ 00. For u = 1, the value 
of Uc is 3.30 which is in good agreement with the value 3.38 found by exact calculation. The 
method is found to work for all values of tu i.e in both the bulk and the adsorbed regimes. 
However, as u is increased the values of Gf,{uJ,u) do not remain as smooth as at lower values 
of uj, therefore introducing some inaccuracy in the values of Uc- The values of {uJc,Uc) at the 
multicritical point found by us are (2.16, 3.30) which are in good agreement with exactly known 
values (2.19,3.38). As 00 and u becomes larger, higher order terms in equation (2.1) become 
significant and therefore to have proper convergence more terms must be considered. This 



explains the difference in ujc2 line found by exact calculation [|T4| and the present calculation 
(see fig. 3). 

III. RESULTS 

The phase diagram shown in Fig. 3 has four phase boundaries instead of three as reported 



in earlier work ||T^,IT5|. The hne Uc separates the expanded and collapsed phases. 



Exact results [14] 
' Our results 




Fig. 3; Phase digram of a directed linear polymer chain in 2-D space. The lo and u axes 
represents, the Boltzmann factor of surface attraction and monomer-monomer attraction respectively. 



Solid line curve corresponds to present calculation and the dotted line to ref. |14|. 



The special adsorption line uj^ separates adsorbed expanded (AE) phase from that of 
desorbed expanded (DE). Beyond 6'-point we have two boundaries uod and ujc.2- The line uod 
separates the desorbed collapsed (DC) bulk phase from that of a surface attached globule 
(SAG) state, whereas the boundary ujc2 separates the SAG phase from the AE phase. The 
point where Uc line meets the special adsorption line Uc-, all four phases AE, DE, DC and 
SAG coexist. The SAG phase which exists between the boundaries Ud and ujc2 for u > Uc 
is essentially a two-dimensional globule sticking to the surface in the same way as a liquid 
droplet may lie on a partially wet surface. These features are in agreement with those of 
isotropic (non directed) polymer chain. The average fraction of monomer on the surface and 
the average number of pairs are calculated using the relations. 

dG , , , aC , 

< Hs >= lim w—— L, and < n„ >= lim u—— L 



The results are shown in Figs. 4 and 5. The transition points on each curve is marked by 
dots. The dot on a curve oi < Ug > (in Fig. 4) indicates adsorption or polymer chain getting 
attached to the surface {i.e u^ and ujd lines depending on the values of -u), the dot on a curve 
< Up > (in Fig. 5) indicates the collapsed-expanded states transition {i. e. point on line Uc 
and iJc2 depending on the values of ou). We may note that the value of < n^ > for SAG phase 
is comparable to that in the DC phase. 
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Fig. 4; Average fraction of monomers on the surface < n^ > as a function of w. Dot on the curve 
represents transition point. 



We have found that along the line Uc and Ud, < Ug >= 0.04 ± 0.01 and along the line u^ 
and LJc2, < Up >= 0.54 ± 0.06. 
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Fig. 5; Average number of pairs < n^, > as a function of u. Dot on the curve represents the 
transition point. 



IV. CONCLUSIONS 



It is obvious from these results that when the chain gets adsorbed from the expanded 
bulk state {i.e for u < Uc), it acquires a conformation at cj ~ ujc{u) such that a small 
fraction of monomers get attached to the surface and others forms a layer parallel to the 
surface. This affects the free energy of the system. At the transition point the polymer chain 
fluctuates between adsorbed and desorbed conformations. Therefore one may get peak in 
< A^^ > — < Ns >^ at uj = ujc which shows divergence as A^ ^ oo. As cj is increased for 
the same value of u, the fluctuations along the normal to the surface get suppressed and at 
large uj{» ujc) the chain lies on the surface with very little fluctuations. On the other hand, 
when the adsorbing chain was in collapsed bulk state, then at cj = Ud the collapsed chain gets 
attached to the surface. Here again the number of monomers getting attached to the surface 
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are small. Since for uo < uod the globule may be any where in the available space and get 
attached to the surface at cj = c^ci, there is a breaking of translational invariance ai uj = ujd 
For ujci < u; < uJc2{u) the chain remains in the form of a globule attached to the surface. 
In this range the monomer- monomer attraction remains effective in holding the monomers 
in the neighbourhood of each other than the surface-monomer attraction whose tendency is 
to spread the chain on the surface. Since dX uo = uod the globule retains its shape there is 
no change in the bulk free energy. Only the surface contribution to the free energy which 
was zero for uo < Ud starts contributing. Therefore, in this case we may therefore not see 
divergence in the bulk specific heat peak in the N —>■ oo limit. For uj > u!c2{u) the globule 
conformation becomes unstable as surface-monomer attraction becomes more effective than 
the monomer-monomer attraction and therefore the chain spreads over the surface (just like 
a liquid droplet spreads over a wetting surface). 

Summarizing we studied a polymer chain immersed in a poor solvent in the presence of an 
attracting impenetrable wall and obtained the phase boundaries separating different phases of 
the polymer chain from data obtained by exact enumerations. We report a state in between 
the desorbed collapsed and adsorbed phases which has conformation of a compact globule 
sticking to a surface in same way as a droplet may lie on a partially wet surface. A similar 
phase also been found in case of isotropic (non directed) polymer chain in both two and 



three dimensions [|Tl|)|lU • The number of monomers on the surface and the number of nearest 



neighbors in different regimes of the phase diagram are obtained. 
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